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Multiple-streaming and the Probability Distribution of Density in 

Redshift Space 

Lam Hui 1 , Lev Kofman 2 and Sergei F. Shandarin 3 
ABSTRACT 

We examine several aspects of redshift distortions by expressing the redshift- 
space density in terms of the eigenvalues and orientation of the local Lagrangian 
deformation tensor. We explore the importance of multiple-streaming using the 
Zel'dovich approximation (ZA), and compute the average number of streams in 
both real and redshift-space. It is found that multiple-streaming can be signifi- 
cant in redshift-space but negligible in real-space, even at moderate values of the 
linear fluctuation amplitude (er^<l). Moreover, unlike their real-space counter- 
parts, redshift-space multiple-streams can flow past each other with minimal 
interactions. Such nonlinear redshift-space effects, which are physically distinct 
from the fingers-of-God due to small-scale virialized motions, might in part ex- 
plain the well-known departure of redshift distortions from the classic linear pre- 
diction by Kaiser (1987), even at relatively large scales where the corresponding 
density field in real-space is well described by linear perturbation theory. We also 
compute using the ZA the probability distribution function (PDF) of density, as 
well as S3, in real and redshift-space, and compare it with the PDF measured 
from N-body simulations. The role of caustics in defining the character of the 
high density tail is examined. It is found that (non-Lagrangian) smoothing, due 
to both finite resolution or discreteness and small-scale velocity dispersions, is 
very effective in erasing caustic structures, unless the initial power spectrum is 
sufficiently truncated. 

Subject headings: cosmology: theory — gravitation — large-scale structure of 
universe 

1 NASA/Fcrmilab Astrophysics Center, Fermi National Accelerator Laboratory, Batavia, IL 60510; e-mail: 
lhui@fnal.gov 

2 CITA, 60 St. George St., Toronto, M5S 3H8, Canada; e-mail: kojman@cita.utor onto. ca 

3 Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045; e-mail: 
sergei@kusmos.phsx.ukans.edu 



-2- 



1. Introduction 

The distortion of the density field in redshift-space by peculiar motion is an old subject 
(see e.g. [Sargent fc 'turner 1977| ; [Bean et al. 1983| ; pavis fc Peebles 1983| ; [Kaiser 1987| ; 



[Hamilton 1992 ). It is commonly held that on small scales, virialized motion causes the 



stretching of structures along the line of sight, creating what is often called the fingers- 
of-God, while on large scales, coherent infall into high density regions causes compression, 
thereby making structures appear more wall-like in the tangential directions than they are 
in real-space. The reader is referred to Hamilton ( |1996|) for a comprehensive review. 

It has also been recognized for some time, however, that the classic linear prediction 
by Kaiser ( |1987| ) overpredicts the compression, even at relatively large scales where the 
corresponding density field in real-space is quite linear (ag^zl). Conversely, the dilution 
effect embodied by the finger s-of- God appears to extend to rather large scales. (See e.g. 
Suto fc Suginohara 1991) ; [Fisher et al. 1993| ; |Gramann et al. 1993| ; Prainerd et al. 1996 



Bromley et al. 1997] .) We will refer to this phenomenon as translinear redshift distortions. 

This effect has often been discussed in the context of attempts to measure (3 (Qf^/b, 
where b is the bias) using the ratio of the real-space to redshift-space two-point function 



or the quadrupole-to-monopole ratio of the redshift-space power spectrum (see e.g. |Cole 



et al. 19941 ). It is found, both from observations, as well as from N-body simulations, that 



the linear theory predictions for these ratios fail even at very large scales, where perturbation 
theory is known to work well in real-space. In general, one finds that the amount of line- 
of-sight squashing of the redshift-space two-point function is overpredicted by the linear 
perturbative formulae. The large uncertainties in published values of f3 are at least in part 
due to our poor understanding of this effect (see Table 1 of Hamilton 1996| ). 



A related phenomenon has also been observed in studies of one-point statistics of the 
density field in redshift-space, albeit in a slightly more subtle form. Hivon et al. (|1995|) com- 



puted the leading perturbative contribution to S3 and S 4 in redshift-space. They observed 
that the agreement on large scales between the predicted values and the measured values 
from N-body simulations is not as good as it is in the case of real-space. Only at very small ae 
(the linear fluctuation amplitude in real-space; in other words, very large smoothing scales) 
did they obtain rough agreements. In fact, above some moderate ae (~ 0.2), the value of 
in redshift-space generally lies below that in real-space, suggesting an effect that resembles 
the finger s-of- God. 



Recently, Taylor & Hamilton ( 1996 ) and Fisher & Nusser ( |1996j ) (see also Hatton & 



Cole |1997| for detailed comparisons with simulations) revisited this issue by computing the 
quadrupole-to-monopole ratio of the redshift-space power spectrum using the Zel'dovich 
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approximation ( Zerdovich 1970 ; henceforth as ZA). Both pairs of authors were able to 



obtain accurate predictions for the shape of the quadrupole-to-monopole ratio on translinear 
scales. In particular, they demonstrated analytically that the departure from the linear 
prediction could take place at relatively large scales, and they interpreted the success of the 
ZA in accounting for this effect as an indication that the failure of the linear prediction on 
translinear scales have more to do with coherent infall rather than virialized motion. 

It is our aim here to continue this line of investigation of nonlinear effects in redshift 
distortions, but we will focus our attention on one-point statistics. Following earlier work by 
McGill ( |1990| ), we explore the role of caustics and multiple-streaming in redshift-space. It was 
pointed out by McGill ( |1990| ), who used one- dimensional Zel'dovich dynamics, that redshift- 
space caustics could form even when the real-space density field is only mildly nonlinear. 
We will quantify this by computing the average number of streams in both real and redshift- 
space using the ZA. As will be shown, at ~ 1, the degree of multiple-streaming could be 



significant in redshift-space but negligible in real-space. This is covered in §3.2 



The above finding indicates that the real to redshift-space mapping enhances the level 
of nonlinearity. An interesting feature of multiple-streaming in redshift-space is that, un- 
like their real-space counterparts, redshift-space multiple streams to a first approximation 
can flow freely past each other (at least before real-space caustics actually form), because 
they arise from physically distinct regions in real-space. This might provide an interesting 
explanation for why the famous Kaiser formulae overpredicts the amount of redshift-space 
compression along the line of sight, even at moderate cr/s. Redshift-space multiple-streams 
could then provide a mechanism for counteracting the linear compression. The free-crossing 
of redshift-streams in the intermediate regime is not unlike the familiar phenomenon of the 
thickening of the ZA pancake. In this sense, the ZA, which can be viewed as a local ap- 
proximation where interaction is ignored (see |Shandarin fc Zel'dovich 1989| for a review; see 



also |Kofman fc Pogosyan 1995| ; [Hui fc Bertschinger 1996| ), is well-suited for studying this 



phenomenon. 

The formalism for computing the average number of streams was in fact originally 
developed to calculate the full one-point probability distribution function (PDF hereafter) 
of density. Our work here is a natural extension of earlier work by Kofman et al. (|1994| ) who 
computed the PDF in real-space using the ZA. We will also use the same methodology to 
compute the PDF as well as the skewness S3 in both real and redshift-space in the small ae 
limit. The formalism is laid out in §|2], the computational method is explained in § |3.1| , and 
the results of the PDF and S3 calculation are discussed in §1373 and 



We particularly focus our attention on the high density regime of the PDF. Catastrophe 
theory tells us that caustics generally induce a p -3 tail, which, however, is not observed 
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in N-body simulations with CDM-type (Cold-Dark-Matter) power spectra, either in real 
or redshift-space. In fact, an intriguing property of the PDF in CDM models is that it 
is well-fit by a single parameter lognormal distribution (in both real and redshift-space), 
where the caustic-tail is absent (for real space, see e.g. [Kofman et al. 1994 ; Bernardcau| 
Kofman 1995[ ). We explore the reason for it in §§, and conclude that smoothing in 
redshift-space (or in real-space), due to both finite resolution or discreteness and small-scale 
velocity dispersions, plays an important role in erasing the caustic-feature in the PDF 0. To 
understand properly the behavior of the PDF then requires a non-local calculation, which is 
beyond the scope of this paper. We briefly discuss possible lines of attack in §[5|. 



2. Formalism 

For the purpose of this article, we denote by p and p z the densities in real (or Eulerian) 
and redshift space respectively, normalized by the mean density (i.e. (p) = (p z ) = 1). The 
evolution of p is determined by the mapping from Lagrangian to Eulerian space, according 
to mass conservation: 



P 



det" 1 



where x l and g J are the Eulerian and Lagrangian coordinates. They are related by the 
displacement 

x = q + D. (2) 



Similarly, the evolution of the redshift-space density is determined by the mapping from 
Eulerian to redshift-space: 

1 + z , 

(3) 



H 



v ■ e e 



where s is the redshift-space coordinate, z is the mean redshift of interest, H is the Hubble 
parameter at redshift z, v is the peculiar velocity and e is the unit vector along the line of 
sight. 



4 Smearing of caustics in the HDM model due to the thermal velocities was considered by Zel'dovich & 
Shandarin Jl982| ) and Kotok & Shandarin ( |1987| ). 
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The redshift-space density p z is then given by (] 



det" 



dqi 



det" 



ds i 



1 + 



1 + z dv 
H dx i 



ee' 



(4) 



where eq. (0) has been used for the second equality We have used the remote-observer 
approximation for the third equality, in other words, setting the term de/dx l = 0. This is 
equivalent to assuming the directional unit vector does not change appreciably for a small 
displacement in the transverse direction. In addition, an interesting matrix identity has also 
been used: det[5y + A^Bj] = 1 + A ■ B, where A and B are arbitrary vectors. 

Let us now introduce the Zel'dovich approximation (ZA) to model the dynamics. 

The ZA consists of assuming D oc D + where D + is the linear time-dependent growth 
factor. Growing mode (in other words, irrotational) initial condition then implies 

B = -D + V^, (5) 

where ip is the displacement potential and V q in component form is d/dq l . 

The peculiar velocity is then simply v — D. where the derivative is with respect to 
conformal time for a fixed q. 



Putting the above into eq. 

Pz = 

where 



Xii — $ii + D 



we arrive at the following expression for p z 
det _1 X (1 + fXr h 1 D k j e i e j )- 1 



D, 



(6) 



(7) 



dq l dqi 

and / = Da/ (Da) ~ (Peebles 1980|) ,p| with a being the scale factor and fl m being the 
fraction of critical density in matter. 

Therefore, for a given Q m , p z is uniquely determined by the line-of-sight unit vector e l 
and the displacement matrix Dij. The probability distribution of p z is then dictated by the 
probability distributions of the later two quantities. 

For random orientations of the line-of-sight, the distribution of e l is given by 

P(e i )de i = sin^, 



1 This is the single-stream redshift-space density. See discussion at the end of the section. 



b A more accurate expression for / can be found in Lahav et al. 1991 and Carroll et al. 1992, but we will 
adopt / ~ ^m 6 i which is a very good approximation at z = 0. 
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where [e x ,e 2 , e 3 ] = [sin#cos0, sin#sin0, cos6*] . 

For Gaussian random initial conditions, the probability distribution, in Lagrangian 
space, of the displacement matrix can be computed exactly. One can always go to 
a frame in which Dij is diagonal, keeping in mind that the line-of-sight is randomly oriented 
with respect to the eigenvectors of this matrix. All we need then is the probability distribu- 
tion of its eigenvalues. We will denote this probability distribution in Lagrangian space by 
P(\j)d\j, where A/s (j = 1,2,3) are the eigenvalues of —D^, Ai > A2 > A3. 

Accordingly, eq. (|5]) can be rewritten as 

Pz = |(1 - Ax)(l - A 2 )(l - A 3 )(i - /E (^y^) 2 )!" 1 (9) 

or further as 

Pz = |(i + / )(i-A 1 )(l-A 2 )(l-A 3 ) (10) 
-f[(l - A 2 )(l - A 3 )e? + (1 - A0(1 - \,)e\ + (1 - A x )(l - A 2 )e*] I" 1 

Thus, the density in redshift-space can ultimately be expressed in terms of the eigen- 
values and orientation of the local Lagrangian deformation tensor. The statistical properties 
of p z can then be inferred from the well-known statistics of the deformation tensor. Eqs. 

(|TT|) represents a redshift-space generalization of the familiar Zel'dovich formulae for the 
density in real-space. Just like the Zel'dovich result in real space, it sheds light on the main 
features of the local structures in redshift space. The first three factors in Eq. (Q) obviously 
arise from the mapping from the Lagrangian to Eulerian (or real) space and the last factor 
comes from that from real to redshift-space. The latter mapping strongly depends on the 
orientation of the deformation tensor with respect to the line of sight e. 

In Fig. H, we show an example of the density evolution of of a fluid element, with the 
eigenvalues today (D + = 1) of Ai = 1,A 2 = 0.5 and A 3 = —0.5. Two cases are considered: 
when the line of sight is aligned with the axis of fastest collapse (Ai) (short dashed line), 
and when the line of sight is aligned with the eigenvector associated with A 3 (long dashed 
line). The evolution of the real-space density is shown as a solid line. The real-space density 
formally diverges at D + = 1 because of caustic-formation. The redshift-space density tends 
to become nonlinear even earlier. In particular, this redshift-space amplification works most 
efficiently if the axis of the fastest collapse is aligned with the line of sight. What we have 
here can be interpreted as a nonlinear manifestation of the linear effect noticed by Kaiser 
(1987). Note how the redshift-space density decreases after (redshift-) caustic formation (at 
D + = 0.5, for the short dashed line). This corresponds to the free-crossing of redshift-space 
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multiple-streams, which counteracts the Kaiser compression of structures along the line of 
sight. 

It is also important to emphasize that eq. @ only gives us a local picture of the density 
distribution. As in the familiar case of the ZA in real-space, our local Lagrangian treatment 
suggests an anisotropic structure, except that the structure is preferentially transverse to 
the line of sight. However, it should be cautioned that the above equations tell us nothing 
about the global distribution of these local structures, whether they are aligned in a filament 
or a sheet is actually determined by the correlation function of shear (see Pond et al. 1996| ) . 



Moreover, small-scale velocity dispersions are likely to at least smear out the caustics, which 
we will come back to later. The only situation in which we can say something about the 
global distribution of structures using our local Lagrangian description is for a truncated 
initial power spectrum, in other words with a displacement field that is smoothed on large 
scales in Lagrangian space, similar to what is done in the truncated Zel'dovich approximation 
( [Kofman et al. 1992| ; |Coles et al. T55g ). 



Let us now define two different probability distribution functions (PDF) of density, 
following the treatment of Kofman et al. ( |1994| ). 

The Lagrangian PDF of (redshift-space) density is: 

Pt(pz) = J S D (p z - p 2 ')P(A,)dA,P(e i ) C /e i , (11) 

where 5d is the Dirac delta function, p z ' in its argument is a function of Xj and e* as 
expressed in eq. (f|), P(e l )de z is given by eq. (||), and P(Xj)dXj is the probability that a 
given Lagrangian element has the set of eigenvalues A/s falling in the prescribed ranges. In 
other words, Ph{Pz)dp z tells us the fraction of the Lagrangian volume that has the redshift- 
space density p z ± dp z /2. 

To obtain the redshift-space PDF of the (redshift-space) density, two effects need to be 
taken into account. The redshift-space PDF is supposed to tell us the fraction of the redshift 
space volume that has the redshift-space density within a particular range, say p z ± dp z /2. 
The first effect arises from the fact that the Lagrangian and redshift-space PDFs differ by a 
multiplicative volume factor: a given Lagrangian volume corresponds to a different redshift- 
space volume. This factor is none other than p z itself. In other words, the redshift-space 
PDF should be equal to the Lagrangian PDF divided by p z . 

If there were no multiple-streaming, this would be the whole story. The second effect 
that has to be taken into account is illustrated schematically in Fig. [|. Multiple-streaming 
is said to occur in places where the mapping from q to s is no longer one-to-one. 

The expression for p z in eq. (f|) is defined for every Lagrangian coordinate q. In 
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other words, it is a single-stream redshift-space density. In regions where multiple-streaming 
occurs, the actual redshift-space density, at a particular s say, should be a sum of p z 's over all 
q's that map onto the same s. Such a non-local calculation is beyond the scope of this paper. 
Instead, following Kofman et al. ( |1994| ), we assume the fraction of the total redshift-space 
volume that is occupied by these multiple-stream-regions is not large (i.e. mildly nonlinear 
regime). However, there are inevitably some such regions. A renormalizing factor N s is then 
needed when computing the redshift-space PDF: 

PApz) = N-'p^P^Pz) , (12) 
where Pl is given in eq. (|TT|) and 

/•CO 

N s = / p x ~ 1 PL(Pz)dp z . (13) 
Jo 

The factor of p z ~ l takes care of the first effect mentioned above: it is the ratio of 
redshift-space volume to Lagrangian-space volume. However, since p z is the single-stream 
redshift-space density, an overcounting of the total redshift-space volume over the total 
Lagrangian-space volume occurs. The true ratio of the two should be unity. N s gives the 
ratio of the total single- stream redshift-space volume to the total Lagrangian-space volume 
(the latter is also equal to the true total redshift-space volume with no overcounting), which 
is in general larger than 1 because of multiple-streaming (or, in other words, overcounting). 
This quantity provides the correct renormalizing factor for P z . 

An additional bonus of the above formalism is that we have a quantitative measure of 
the degree of multiple-streaming in N s . It tells us the average number of streams at an 
arbitrary point in redshift-space. Fig. ^| provides a nice illustration of how N s should be 
interpreted. Suppose V T is the total true (no overcounting) volume in redshift-space (i.e. 
the "volume" shown on the y-axis). Suppose further that V 1 is the part of V T that is in 
single-stream regions, and V 3 is the part that is in triple-stream regions i.e. V T = V 1 + V 3 . 
Recall that N s is the ratio of the total single-stream redshift-space volume to the total true 
redshift-space volume, which means N s = (V 1 + 3V 3 )/V T = P(l) + 3P(3), where P(l) and 
P(3) are the probabilities that a given region is in single-stream and triple-stream regimes 
respectively. Hence, N s is exactly what one would refer to as the average number of streams 
at a point in redshift-space. Obviously, the argument extends to higher (but always odd) 
number of streams. 

Note furthermore that while eq. ([12]) gives only an approximate redshift-space PDF in 
the limit of little multiple-streaming (N s close to 1), eq. (^) is an exact expression for the 
average number of streams, at least within the framework of the ZA. This should be kept in 
mind in the discussion of the next section. 
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Lastly, note that the above expressions (eq. 0, [|TTJ, [|12|], reduce to their real-space 
counterparts in the limit / = 0, in other words, when the real to redshift-space mapping is 
trivial. 



3. The Average Number of Streams and the PDF: Method and Results 

3.1. Method 



The quantities we are interested in are P z (p z ) arid N s in eq. (|TJ) and (|HJ). Both of them 
require performing the integral in eq. (fTTl). Kofman et al. ( |1994j ) succeeded in reducing the 
analog of this integral in real-space (i.e. setting / = in eq. 0) to a one-dimensional integral 
which has to be done numerically. Redshift-space distortions introduce extra complications. 
We will instead perform the integral in eq. (jTTf) numerically right from the beginning. The 
delta function restriction 5d[Pz ~ Pz'(^j,£i)] is most easily dealt with using Monte Carlo 
techniques (see Keiscncggcr fc Miralda-Escudc 19951 f° r a similar calculation in the context 
of the PDF of the optical depth rather than the density in redshift-space). In other words, 
we generate realizations of Aj-'s and e*'s, and identify combinations of them which satisfy 
the delta function restriction. The directional unit vector e l is easily handled because it is 
randomly distributed over the 2-sphere, with the probability distribution given by eq. @. 
On the other hand, the probability distribution for the A^-'s is somewhat more complicated. 

The probability distribution of the eigenvalues of the displacement matrix Dij, for a 
Gaussian random field, was first calculated by Doroshkevich (|1970|) . We adopt a form for 
this distribution which is better suited for the Monte Carlo computation outlined above 
(following |Bertschinger fc Jain 199~4] ; see also [Reisenegger fc Miralda-Escude 1995| ). Let 
— Aj's be the eigenvalues of the displacement matrix D^, and let A, 7 and a be variables 
defined through A/s as follows 

A A a + 27T a — 2n a 

[Ai, A 2 , A 3 J = — cos7[l, 1, 1J + -^=sni7[cos( — - — ), cos( — - — ), cos(-)J , (14) 

where the A ranges from to infinity, and both 7 and a are between and it. The above 
ranges enforce the ordering Ai < A2 < A3. It can be shown that for any given combination 
of Aj's in this order, there corresponds a unique set of A, 7 and a. 

The probability distribution of X/s is given by ( [Bertschinger fc Jain 199^ ): 



2 r 



n 2\ 



2<7J 



dA 2 



— sin 4 7(i7 

37T 



-sinada 
2 



(15) 



where 01 is the root-mean-squared (rms) amplitude of the linear real-space density fluctua- 
tion. 
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A nice feature of the set of variables A, 7 and a is that their joint probability distribution 
factorizes. A 2 is distributed like a gamma-deviate. 8 sin 4 7/(37r)c?7 can be integrated exactly 
to give drj/2 where 77 = 1 — 27/7? + 4sin27/(37r) — sin47/(67r), so that 7] is distributed 
uniformly between —1 and 1. Lastly, cosa is also uniformly distributed between —1 and 1. 
Standard numerical methods exist for generating variables distributed in the above manner 
(see Press et al. 1992]) . 

Note that the only free parameters in the above calculation are <jg and Q m , or more 
precisely, /. The former controls the distribution of X/s through eq. (|15|), while the latter 
determines the degree of redshift-space distortion (eq. or [[|). 



3.2. The Average Number of (Lagrangian) Streams 

Fig. ^ shows the result of the calculation outlined above for N s , the average number of 
Lagrangian streams, defined in eq. (|13|). This is shown as a function of ae (eq. |IIJ), the rms 
amplitude of the linear real-space density fluctuation, for three different values of fl m . The 
Q m = case corresponds to the average number of streams in real-space. 7 

Several points should be noted. First, redshift-space distortion, or peculiar motion, has 
the effect of increasing the amount of multiple-streaming over that in real-space. Second, a 
larger Q m means more redshift-space distortion, hence a higher average number of streams. 

At oi = 1, the average number of streams in real-space is very close to 1.0, while 
the average number of streams in redshift-space is 1.1 and 1.2 for fl m = 0.3 and fl m = 1 
respectively. This means multiple-streaming can be insignificant in real-space, but going on 
in earnest in redshift-space.0 In other words, even at relatively small ae, when the real-space 
density field is quite linear, redshift-space distortion introduces inherently non-perturbative 
effects through multiple-streaming. 

The above result perhaps explains the success of perturbation theory, which is by 
nature a single-stream theory, in estimating various statistics of the density field in real- 



space, at op^l (Many references should be cited here, see LJuszkiewicz fc Bouchet 1995 



and Pernardeau 1996| for short reviews. See also [Scoccimarro fc Frieman 1996| and ^coccT 



^Fig. 1 of Kofman et al. ( 1994 ) contains an error for the real-space average number of streams. For the 



perturbative result in real-space, see also Munshi, Sahni & Starobinsky (1994). 

8 1.2 might not seem too different from 1.0. However, it should be kept in mind this means, at any 
given point in redshift-space, one finds on the average 1.2 streams. In other words, the total single-stream 
redshift-volume is 20% higher than the total true redshift- volume. 
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marro 1997 for loop corrections which can be important at large enough <?£.). The same does 
not appear to be true for statistics of the redshift-space density field, however. For instance, 
it is well-known that the linear-theory redshift distortion factor deduced by Kaiser ( J1987 ) 



fails to describe the quadrupole-to-monopole ratio of the redshift-space power spectrum even 
at very large scales (see Hatton & Cole |1997| and references therein for a recent discussion). 



This failure is an indication of the importance of non-perturbative effects such as those due 
to multiple-streaming in redshift-space. 

The fact that multiple-streaming in real-space is less severe than multiple-streaming in 
redshift-space has another interesting implication. Recall the Lagrangian to real-space to 
redshift-space mapping: q — > x — > s. The multiple- (Lagrangian)-streaming in redshift-space 
arises not so much from the q — > x mapping (otherwise, the average number of Lagrangian 
streams in real-space would not be so close to 1 for, say, og = 1), but rather from the x — > s 
mapping. 

This is interesting because the multiple-streams at a given s then mostly come from 
distinct places in real space i.e. they do not interact strongly, unlike the case of multiple- 
(Lagrangian)-streams at the same x. One well-known failure of the Zel'dovich Approximation 
is that it does not account for the gravitational interaction between different Lagrangian 
streams around a real-space caustic i.e. the different streams (unphysically) freely go past 
each other at a Zel'dovich pancake. Around a redshift-space caustic, however, the multiple 
streams can (approximately) free-stream past each other because they are not located at 
the same real-space location. This might provide a mechanism for partially canceling the 
line-of-sight compression of structures pointed out by Kaiser (1987). 



3.3. The PDF of Redshift-space and Real-Space Density 



We compute the PDF of redshift-space density according to eq. (|12|) for Vt m = 1 and 
Q m = 0, and for 2 different cr^'s: the result is shown in Fig. f|. The Q m = curves correspond 
also to the PDFs of the real-space density. We have chosen a/s small enough so that the 
average number of streams is close to 1 (see Fig. |^). This is to ensure the validity of the 
approximation involved in eq. flL?]). 



It can be seen that peculiar motion has the effect of making the PDF look more nonlinear 
in redshift-space than it is in real-space. This should be expected based on the Kaiser 
( Kaiser 1987|) effect alone: large scale coherent infall into initial density enhancements (in 
real-space) shifts such regions further out into the high-density tail of the PDF, while at the 
same time creating lower-density voids thereby making the PDF peak at lower densities. 
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At higher cr/s, translinear distortions, or redshift-space multiple-streaming, would start 
to cancel the Kaiser compression, and begin to cause dilution around the initial density 
enhancements. At still higher cr/s, a qualitatively different dilution effect, the finger-of-God, 
kicks in due to the formation of virialized clusters. This happens after multiple-streaming 
becomes important in real-space. Unfortunately, the PDF as expressed in ([|) is incapable of 
taking multiple-streaming (in either redshift or real-space) properly into account, as we have 
explained in §|2|. Nonetheless, we will attempt to probe into the regime of multiple-streaming 
by considering the effect of caustics on the high p z tail of the PDF in §|4]. 

How good is the ZA-predicted PDF as a description of the true PDF? We show in Fig. 
|5] a comparison of the ZA PDF versus the PDF measured from an Q m = 1 CDM simulation, 
for two different output times, smoothed at the same comoving scale. The linear CDM power 
spectrum has T = 0.25, and the box size is 100 h _1 Mpc, simulated on a 128 3 grid. The reader 



is referred to the Hydra consortium flhttp : // coho . astro . uwo . ca/pub/ data . html] ; see also 



Couchman et al. 1995|) for further details of the simulation. The one free parameter we have 



to fix for the ZA prediction is the linear fluctuation amplitude oi. It can be seen that, as 
expected, ag is always smaller than the fluctuation amplitude of the evolved-and-smoothed 
density field, which we call a Q \. We also show lognormal distributions which reproduce the 
evolved fluctuation amplitude of the simulation outputs. In other words, 

Pio g n.(p z ) = -f= e , (16) 



where cr e g is chosen to equal yln(l + cr^). It is impressive how well it fits the N-body 
data out to large p 2 's. The same phenomenon has been observed in real-space for a similar 
cosmological model ( [Kofman et al. 1994| ). The error-bars are estimated by dividing the box 



into 4 different subsamples, and computing the dispersion. It should be emphasized that we 
have purposefully chosen outputs, with suitable smoothing, that are only mildly nonlinear, 
so that the average number of Lagrangian streams is very close to 1, making eq. flT2"| ) a good 
approximation. 

The overall agreement between the ZA PDF and the N-body data is decent at these 
low values of a?. However, the ZA PDF tends to overpredict the probability of high p z . It is 
in fact well-known that the ZA PDF has a long tail of P{p z ) oc I /pi, due to the presence of 
caustics. This is true in both real and redshift-space. It implies that, unfortunately, oi has 
to be left as a free parameter in the above comparison because the ZA PDF does not yield 
a well-defined cr nl . We will revisit this issue in §0L 
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3.4. S 3 in Redshift-space and Real-space 



s. = < 17 > 



The methodology outlined in § |3.1| can be easily adopted to compute other one-point 
statistics, such as S 3 : 

{{P Z - I) 3 ) 
<(P. " I) 2 ) 

As in the case of cr nl , because of the caustic-induced tail of the ZA PDF, S 3 is, strictly 
speaking, undefined except in the limit of vanishing <j(_. The finite number of realizations 
used in our Monte Carlo method provides a natural regularization by truncating at some 
high p z 's, thereby giving a finite 5*3. As one decreases the input parameter ag in the Monte 
Carlo integration, the output S 3 converges eventually to the correct ZA predicted value 
for infinitesimal ag. Our results are summarized in Table [I] (see also e.g. Pernardeau fc 



Kofman 1995| for the real-space ZA prediction). We find that the rate of convergence for S4, 



S5, etc becomes progressively slower, and so we only give the 5*3 values here. 

It can be seen that peculiar motion tends to slightly increase S3 over its real-space value, 
as the reader might have guessed from the PDFs shown in Fig. |4j. This behavior is consistent 
with what is found by Hivon et el. ( |1995| ) in the case of the exact dynamics, in the small og 
limit. However, as is well-known in the case of real-space, in the perturbative limit, the ZA 
systematically gives lower S3 than the exact dynamics. 



As emphasized by Hivon et al. ( |1995[ ), the perturbative prediction breaks down even 



at relatively small c^'s (say 0.2), or large smoothing scales, and the actual N-body S3 is 
generally smaller in redshift-space than in real-space. They attributed it to some kind of 
finger-of-God effect on large scales, which we view as in part due to multiple-streaming in 
redshift-space (that is distinct from the familiar small-scale finger-of-God). Unfortunately, 
the ZA yields a meaningful skewness only in the vanishing <jg limit, making it difficult to 
investigate the impact of these effects on S3. 



4. Erasing Caustics ? 

Why does the ZA PDF predict a long high-p^ tail, and why is it not observed in N-body 
simulations, at least for CDM-type models? We believe the answer to the latter at least in 
part has to do with smoothing due to finite resolution or discreteness as well as small-scale 
velocity dispersions. But let us first briefly review where the high p z tail comes from. The 
reader is referred to Shandarin & Zel'dovich ( |1989| ) for a review. Most of our arguments 
below apply equally well to real-space caustics. 

A caustic is located at a point where p z or det _1 9sj/ 'dqj diverges. As Zel'dovich (|1970| ) 
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originally argued, this is generically due to the vanishing of only one of the eigenvalues of 
the matrix dsi/dqj. Hence, we can focus our attention on the direction that aligns with 
the eigenvector associated with the vanishing eigenvalue. We will denote the relevant one- 
dimensional coordinates simply by s and q. The situation is depicted in Fig. |2|. 

Taylor expanding s as a function of q around q = q c (corresponding to s = s c where the 
caustic is), one can see that to the lowest order, s is a quadratic function of q i.e. 



Sc ~2dq^ 



q=q c 



In other words, p z should fall off as l/(q — q c ) or l/y/s — s c , around the caustic. Note 
that, strictly speaking, p z should be given by the sum over the contribution from each 
multiple-stream around the caustic (as in the case depicted in Fig. |2|), but the singularity 
due to l/(q — q c ) or 1/y/s — s c dominates. Note also that this singular behavior of p z occurs 
only on one side of the caustic in s-space. 

To understand the caustic's influence on the PDF, the crucial point to remember is that 
one has equal probability of locating at any position in s-space with respect to the location 
of the caustic. Therefore, P(s)ds oc P(p z )dp z , with P(s) independent of s, from which the 
P(Pz) oc 1/pl behavior follows. The premise is that asymptotically high p 2 's occur only 
around caustics. 

Why is the caustic-induced feature of the PDF then not observed? 

We believe the answer lies in final (as opposed to initial or Lagrangian) smoothing. 
There are two sources of final smoothing: one due to finite resolution or the discrete nature 
of the simulations (or observations), and the other due to smearing by small-scale velocity 
dispersions. Two conditions have to be met for the caustic-induced high p z tail to survive 
these kinds of inevitable smoothing in s-space. 

First, it is important that the scale of the caustic in s-space, As caus ti C , is larger than 
the smoothing scale, As smoot h, where the smoothing scale can be taken to be whichever is 
dominant of the two kinds of smoothing mentioned above As smoot h = max(As res ., a v (l + 
z)/H). The smoothing scale As res . is associated with the finite final resolution, and a v is the 
velocity dispersion along the line of sight. 

The scale of the caustic in s-space is related to its counterpart in g-space by eq. (|T8|): 

1 



As caus tic 2 



d 2 s 



dq 



where Ag caust i C is the scale of the caustic in g-space. A good definition for Ag caust i C is given 
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by the magnitude of q — q c such that the next higher order term ignored in eq. becomes 
significant. In other words, Ag caustic ~ [d 3 s / '<9g 3 ] _1 [<9 2 s / 'dq 2 }. One can use even higher order 
terms for its definition, but it would not affect our arguments below. 

Let us then express the first condition for the survival of the caustic-induced feature as 



As sm0 oth *C ^ 



d 2 s 



dq 2 



(20) 



<Jc 



The second condition has to do with the density profile of a smoothed caustic. As- 
sume that the first condition above is met, the smoothed profile of the caustic would 
be given by p z (s) = f ( s ) / \J \ d 2 s / dq 2 \ qc A s smoo th where f(s) is some dimensionless profile 
which, at its peak, is equal to 0(1). This is a very generic prediction quite indepen- 
dent of the form or origin of the smoothing. This means the smoothed caustic peaks at 



Pzmax ~ ^/y\9 2 s/dq 2 \ qc A Ssmooth- Now, if one would like such a smoothed caustic to con- 
tribute to the high p z tail, it is important that p Zmax ^> 1. Hence, the second condition for 
the survival of the caustic-induced feature in the PDF is: 



d 2 s 



dq 2 



1 



As sm0 oth 



(21) 



The two conditions in eq. (|20|) and (^) together imply the necessary requirement that 

As sm0 oth 



Aq, 



< 1. 



(22) 



caustic 



Therefore, the only remaining question is what the typical Ag caustic should be. This 
can be estimated using the definition suggested before: that Ag caust i C is given by [d 3 s / dq 3 ]^ 1 
[d 2 s/dq 2 ] evaluated at the caustic. One can calculate exactly what the average value of 
this ratio should be, for Gaussian initial conditions, imposing the constraint that the first 
derivative vanishes. All one needs is the joint probability distribution of all the derivatives 
up to the third one at a single point. However, there is no need to do any calculation. The 
only relevant scale in this problem is the initial smoothing scale of the density field (or the 
displacement field) in g-space, let us call it Ag smoo th- Therefore Ag caust i C ~ Ag smoo th- 

For CDM-type power spectra, the initial smoothing scale Ag smoo th is of the order of a 
few grid-spacings (unless the initial power spectrum is deliberately smoothed or truncated), 
while the final smoothing scale As smoo th is at least also a few grid-spacings (and probably 
even larger due to small-scale velocity dispersions). This implies the condition in eq. ( |22"D 
is never met, hence explaining why the theoretical caustic-induced tail of the PDF is not 
observed, at least in CDM simulations. In other words, to understand the high-p z behavior 
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of the PDF properly, both multiple-streaming and smoothing have to be properly taken into 
account. This also points to a weakness of the calculation presented in §fj. Because final 
smoothing in s-space is not taken into account in the formalism, effects like what is discussed 
above are missed. We hope to pursue improvement of the calculation in this direction in the 
future. 



5. Discussion 

Using the Zel'dovich approximation, we find that peculiar motion significantly raises 
the degree of multiple-streaming. We quantify this by computing the average number of 
Lagrangian-streams in both real and redshift-space. As is shown in Fig. |||, multiple- 
streaming can be insignificant in real-space, but non-negligible in redshift-space, even at 
moderate values of a^. This implies that most of the multiple-streaming in redshift-space 
arises from the x — > s mapping (real or Eulerian to redshift-space), rather than from the 
q — > x mapping (Lagrangian to Eulerian space). 

That the level of nonlinearity is amplified by the real to redshift-space mapping implies 
non-perturbative effects have to be taken into account in the calculation of redshift distor- 
tions, even on scales where the real-space density field is relatively linear. An interesting 
feature of redshift-space multiple-streaming is that redshift-space multiple-streams can freely 
flow past each other, unlike their real-space counterparts. This counteracts the familiar large 
scale compression of structures along the redshift direction first pointed out by Kaiser (1987), 
but it is physically distinct from the stretching of structures due to virialized motions on 
small scales. This might offer an explanation of why, even on relatively large scales, linear 
theory overpredicts the amount of reshift-space squashing of the two-point function along 
the line of sight. To make this quantitative, however, requires a detailed calculation, which 
we will discuss briefly below. 

We have computed also using the ZA the probability distribution function (PDF) of 
density, as well as S3, in both real and redshift-space. At the largest scales, or small 01 
limit, S3 does not appear to change significantly from real to redshift-space, consistent with 
the finding of Hivon et al. ( |1995| ). We have compared the ZA PDF with that measured 
from a CDM N-body simulation (Fig ^). We find that the PDF in redshift-space is well- fit 
by a lognormal distribution, similar to its counter-part in real-space (|Kofman et al. 1994; 



Bernardeau fc Kofman 1995|) . This suggests the PDF in both real and redshift-space can be 
described by a single-parameter family of functions (<7 e ff in eq. [[R|). This is certainly worth 
exploring with more N-body simulations and different cosmological models (see [Scherrer~%; 
|Gaztanaga 1998| ). 



-17- 



The simulation PDF reveals no caustic-induced high p z tail, contrary to the ZA pre- 
diction. The same has been observed before in real-space as well ( |Kofman et al. 1994| ). We 



argue in §^ that this is at least in part due to the fact that final smoothing in redshift-space 
(or real-space) is not taken into account in our formalism. Smoothing naturally arises be- 
cause of finite resolution or discreteness, as well as due to small-scale velocity dispersions. 
We give the criterion for the survival of the caustic feature through smoothing in eq. (f22| ) : 
the final-smoothing scale has to be smaller than the Lagrangian-smoothing scale. 

In practice, we expect caustics to be efficiently erased by final-smoothing, unless the 
initial power spectrum is sufficiently smoothed or truncated. This is reminiscent of a related 
situation in gravitational lensing: that caustic features can only be observed if the source-size 
is sufficiently small. 

As mentioned before, a proper treatment of the statistics of the density field in redshift- 
space requires tackling the issues of multiple-streaming and final-smoothing. While such a 
calculation is beyond the scope of the present paper, Taylor & Hamilton ( |1996| ) and Fisher 
& Nusser ( |1996| ) have already made significant progress in using the ZA to compute the 
quadrupole-to-monopole ratio of the redshift-space power spectrum. They use the following 
expression for the Fourier transform of the redshift-space density: 

p z (k) = [ rf 3 ge- lk (q+Ds) (23) 



where D s as a function of q is the displacement from the Lagrangian position q to the 
redsfhit-space position s. It comes from the following basic equation for the redshift-space 
density: p 2 (s) = / d 3 g5 j ci(s — q — D s ), which implicitly allows for multiple-streaming by 
automatically summing over density contributions from all possible streams. The key here 
is that the mapping from real to redshift-space is treated exactly, instead of perturbatively, 
hence the non-perturbative nature of the redshift-space density field is taken into accout even 
if the real-space density field is still quite linear. These authors find that the ZA predicts 
the right-shape for the quadrupole to monopole ratio but underestimates the zero-crossing 
scale (see also Hatton & Cole [19971) . One can take their treatment one step further by using 
approximations that are closer to the exact dynamics (see e.g. Scoccimarro |1999~1) . Whatever 
the approach, it appears likely that, because of the prominence of multiple-streaming in 
redshift-space, a pure perturbative calculation would be inadequate. Non-perturbative effects 
such as those due to multiple-streaming (or even pre-multiple-streaming) have to be taken 
into account explicitly. 

Lastly, our Lagrangian analysis here still leaves open the question of how the local 
pancakes or caustics are spatially distributed globally, whether they preferentially lie along 
filaments or sheets. Calculations have shown that filamentary structures in real-space are 
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preferred for Gaussian initial conditions (e.g. Bond et al. 1996 , Yess fc Shandarin 1996| ). 



It remains to be seen how the redshift-space mapping might alter this picture. Our local 
calculation, however, does provide a valid description of the global distribution of structures 
if the initial power spectrum is sufficiently smoothed, as is done in the truncated Zel'dovich 
approximation. 
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S 3 


0.01 





4.0 


0.01 


0.3 


4.1 


0.01 


1.0 


4.2 



Table 1: The ZA predicted S3 in real-space (the Vt m = case) and redshift-space, for 
a t = 0.01 




Fig. 1. — The density evolution of a fluid element, as a function of the growth factor D + . 
The eigenvalues of the deformation tensor (at D + = 1) are Ai = 1, A 2 = 0.5 and A 3 = —0.5. 
The solid line shows the evolution of the real-space density. The short and long dashed 
lines show how the redshift-space density evolves if the line of sight aligns with Ai and A 3 
respectively. 
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Lagrangian space q 



Fig. 2. — A schematic diagram illustrating the mapping from Lagrangian-space to redshift- 
space. Multiple-streaming occurs in places where the mapping is not one-to-one, in other 
words, several q's are mapped to the same s. 




Fig. 3. — The average number of Lagrangian streams in redshift space (N s ) versus the rms 
amplitude of the linear (real-space) density fluctuation (erg), for three different values of Q m . 
The Vt m = curve gives also the average number of Lagrangian streams in real-space. 
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Fig. 4. — The ZA prediction for the PDF of redshift-space (solid line) and real-space (dashed 
line) density for two different amplitudes of the linear (real-space) density fluctuation, Of = 
0.13 and at = 0.49. Q m = 1 for the redshift-space curves. 
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Fig. 5. — Comparison of the ZA-predicted PDFs (solid lines) with the measured PDFs 
(points with error-bars) from CDM, Q m = 1, N-body simulations. The dotted lines represent 
a lognormal fit to the N-body data. The upper panel is from a simulation output at a = 0.23, 
and the lower panel is from an output at a = 0.49, each smoothed with a Gaussian window 
of 6.25 h _1 Mpc. The actual rms amplitudes of the density fluctuation (cr nl ) are 0.34 and 0.49 
respectively, while the linear cr/s used for the ZA-curves are 0.13 and 0.35. 



